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An integral equation for the electric field in a continuously stratified ionized medium 
is derived; this is then manipulated to yield equations for the reflection and transmission 
coefficients, at the same time being susceptible to physical interpretation. The equations 
are solved by successive approximations, the first terms being Fresnel-type coefficients. 
Various applications of the results are discussed. 

1. Introduction 

In a recent paper, Wait [1962] has followed an idea of Brekhovskikh [1960] and has obtained 
approximate forms for the reflection coefficients of VLF and ELF radio waves when the 
ionosphere is not sharply bounded. In the isotropic ease, the method is to obtain a first order 
nonlinear differential equation for l\(z), a function of the height z that reduces to the reflection 
coefficient below the ionized layer. This equation is solved by successive approximations, 
the first approximation being the familiar Fresnel reflection coefficient. Neither Brekhovskikh 
in his comprehensive text (p. 220) nor Wait discusses any physical interpretation of the cor- 
responding integral equation. Here, we show that by means of a different, yet ultimately 
equivalent, formulation of the problem, the resulting equations are susceptible to physical 
interpretation. The theory is applicable to all angles of incidence and to anisotropic iono- 
spheres. The interpretation is made possible by the recent investigations of Westcott [1962a] in 
which he has derived various expressions yielding the effect of particular regions of the 
ionosphere on the reflected wave. 

2. Isotropic Model 

Vertical propagation in an isotropic ionosphere, including the effect of the collision fre- 
quency, is governed by the differential equation 



'M't^* 



dz 2 



where the 2-axis is vertical, E denotes a particular horizontally polarized component of the 
electric field, k=u/c the time factor e io>t being suppressed, X=Ne 2 /e moo 2 , and Z=v/o). We shall 
write l—X/(l—iZ) = 7i 2 , where n is the refractive index. 

We consider the following configuration for the ionized medium. For z<^a there exists 
free space for which n= 1 ; for z^>b there exists homogeneous medium extending to infinity 
for which X=X 1} Z=Z U n = Ui. For a<2<6 the properties of the medium are functions of z, 
varying from n=l to n=Ui. 

Let f denote an intermediate height in the range a<£<b. We multiply the equation 

E"+Jc 2 n 2 E=0 (1) 

by 0~ ,< * w i ( '~f ) and integrate from f to b. To this, we add the result obtained by multiplying 
by e mz ~ v and integrating from a to f. We have, upon integrating by parts twice in succession, 
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fV , «-***t c, -P(fo= J B , ^- tta i 0r -» +ikn 1 f E f e~ ikn ^ z -Vdz 

= (£ ,, +^7i 1 £)6- an i (2 -f ) r~F^f 6 £ , e- aw i (z -f ) ^. (2) 
Ir Jr 

Similarly, 

f r jB? // « tt < , -» ) &=(S / -iikE)« tt(i -f> f -P {*Ee mz ~Vdz. (3) 

We now consider a wave incident below the ionosphere. We therefore restrict ourselves 
to that solution of (1) for which 

E=e-** 9 +Re ik * 0<a), 

E=Te~ ikn i z (z>b), 

z=0 being the phase-reference level for both coefficients R and T. Then 

(E'+ikn^e-^'-v^O, 
and 

(E f -ikE)e mz -^^a=(-ike- ika +ikEe ika -ike- ika -ikR 

=—2ike- ik t. 

The sum of the two integrals (2) and (3) is now given by 

C t E"e mz ~t ) dz+ f E n e~ ikn i (z -^dz 

= -ikE({)-ikn l EU)+2ike- ik !-k 2 \*Ee mz ~Vdz-k 2 n\ f Ee~ ikn ^ z -Vdz. 
Hence, when the final term in eq (1) is taken into account, we obtain 

Ja Jf 

reducing to 

(4) 

This constitutes a suitable integral equation for the field at any point f in the medium; 
the integration process is extended over those portions of the ionized medium where the re- 
fractive index differs from unity or from n 1 . In particular, if free space exists above z—h 
(that is, if 7fci=l), the integral equation becomes 

E(0 = e- iH +lik C Y^ E ^ emz ~^ dz +l ikf Y^EWe-'K'-Vdz. (5) 

It should be pointed out that the right-hand sides of both eqs (4) and (5) are independent 
of a, provided that a is any point in the free space below the ionized region. 

3. Reflection and Transmission Formulae 

In eq (4), let f=a, in which case 
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E({)=e- tka +Re ika . 
Then 

e -ika +Ee i*a=2_ e -i*a + _!* C b fX__X\\ E(z) e - ikn ^~ a) dz 
1+% 1+7?! J a \l—lZ 1— ^Z 1 / 

reducing to 

In particular, if rii=l (that is, if X 1 =0) y this reduces to 

B=i it f y^ E{z)e~ ikz dz. (7) 

It is obvious that the limits may be replaced by — oo and + °o respectively, since the integrand 
vanishes in these extended regions. 

Similarly, we place f = b in eq (4) , in which case 



Hence, 
reducing to 



1+% 



E{{)=Te- ikn i\ 

J % jL- \ ^ F(z)p ik{z ~ b) dz 

+ l+rhJa l-iZ* [Z)e dZ > 



e -im- H )i> + _™ e -«(i-» 1 )» _^ E(z)e ikz dz. (8) 

1 +7&i ,Ja 1 ^^ 



1+% 

If free space exists above z=b, this result simplifies to 

T=i+- 2 ik £y^2 E ^ eikZ(lz - (°) 

It should be pointed out that results (6), (7), (8), and (9) are independent of a and b provided 
that a and b lie in the respective homogeneous regions below and above the ionized medium. 

4. Varying Region Extending to Infinity 

If there exists no homogeneous medium of refractive index Ui above z=b, then the varying 
medium will extend to infinity. In particular, the medium becomes infinitely overdense 
there if X-^co as z->co . Under these circumstances, we multiply eq (1) by e ik{z ~^ and in- 
tegrate from a to f as before, but now we multiply eq (1) b}^ e~ ik{z ~^ and integrate from f to 
oo . The same analysis as before leads to the consideration of 

(#'+*£>" ^"^ I 00 . 

It should be observed that this factor does not converge at the upper limit if a homogeneous 
medium extends to infinity. It was for this reason that previously we employed the factor 
e~ ikn i (z ~V to secure a result that vanished when z = b. For convergence in the present case, 
we require E'+ikE-^O as z-^co . This limit certainly exists if the region becomes infinitely 
overdense, since the allowed W.K.B.J. solution for E then becomes exponentially small in 
magnitude as z— >co . We obtain the result 
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We could have multiplied by e~ ikn ^ z ~^ in the second integral, where n x is arbitrary (Rln{^> 0), 
but this would have been an irrelevant procedure in this case since U\ would have no definite 
meaning in the ionized medium. 
The formula for R becomes 

R=l ik r -X-E(z) e~ ik2 dz, (10) 

where a occurs below the ionization. No formula for T exists in this case. 

5. Physical Interpretation 

There are three distinct contributions to the field E(£) in eq (4). These are: 

Ci)2«-«f/(i+n 1 ); 

(ii) elements of the form — — \^ E{z)8ze mz ~ v originating below the height f; 

(hi) elements of the form— — (- r~— -:: — ^W ) E(z)8ze -**»i<*-f> originating above the 

height f. 

The vector sum of all these waves yields the total electric field at the height f. (i) is a 
kind of transmitted wave complete with the appropriate Fresnel transmission coefficient, 
(ii) represents a plane wave propagated vertically upwards as in free space (according to the 
factor e~ ik *), originating from the elementary layer of thickness 8z situated at 2<£. All such 
layers below f yield waves of similar forms, (hi) represents a plane wave propagated vertically 
downwards as in a homogeneous medium of refractive index n x (according to the factor e ikn ^), 
originating from the elementary layer of thickness 8z situated at z>f. All such layers above 
f yield waves of similar forms. 

The work of Westcott [1962a] shows that this is the correct interpretation of the elements 
involved. Each individual layer of thickness 8z in an ionized medium gives rise to reradiated 
waves, their originating strength depending on the exact value of the electric field within the 
layer 8z. In Westcott 's series of papers he considered these reradiated waves as propagated 
in free space, thereby excluding the possibility of a homogeneous medium existing above z—b 
of refractive index other than unity. His formulae would not converge if such a medium 
existed there. We have modified the theory to allow for this, at the same time providing 
the reflection and transmission formulae derived in section 3. It should be observed that 
our method has not proved the interpretation placed upon the elements involved; to accomplish 
this would require the more elaborate analysis of Westcott, considering from first principles 
the Hertz vectors involved in the reradiation processes. 

6- Approximation Methods 

In order to calculate the reflection coefficient (6), it is necessary to know E{z) throughout 
the medium beforehand. If E(z) were known, R would also be known, implying that there 
would be no need to evaluate R from formula (6), which after all is but an identity. But an 
iterative method of solution maybe adopted, using approximate solutions derived from eq (4). 

A series of approximations may formally be written down based on successive substitution 
into eq (4). We write 



Now 

l#o(r)l=r 



ll+%l 

and if throughout the range a, to b 



max 
then 



Ul-iZ\\\l-iZ \-iZj e |J m > 



M|<^muM«)l 

It follows that the series for E converges rapidly provided 

2kM(b-a) 



1+fli 



-<1, 



that is, provided the layer is sufficiently thin [(b — a) small enough] or provided M is small 
enough (the ionization density small throughout). 

If free space exists above z=b, the formulae simplify. The successive approximations are 

TO =«-**, 

=\ ik£j^ 2; E n - 1 (z)e-'^-i\dz. 

If M denotes the maximum of X/\ 1— iZ\ , the series converges rapidly if 

KfcAf(6-a)<l. 

If X contains a small constant multiplier a, evidently this process yields E(£) as a power-series 
expansion in powers of a, since E n (£)cca n . 

The meaning of the individual terms in the series becomes clear. E^Cf) represents the 
reradiated field if the cause of the reradiation process at each level is given by -Eo(f) and not 
E(£). Similarly, E 2 ({) represents the reradiated field due to the field 2?i(f) only acting on the 
free electrons, and so on. 

Substituting the successive terms of E(£) into eq (6) for li, we obtain the corresponding 
series for R. The first two terms of (6) become 



R^^Ae-^+j^K^e- 






1+th ' (1+n,) 2 

If »i=l above z=b, this reduces to 

^=2 ik £A e ~ 2tkldz (11) 

from eq (7). Under these circumstances, eq (9) yields for the first term in T, 

1 C b X 
T = 1+ 2 ik ) a l=iZ dz - 

This value of R arises from the contributions to the reflected wave produced by the incident 
field e~ ikz being regarded as the sole cause of the reradiation process at each level. It corre- 
sponds, in fact, to the use of the Born approximation. 
For this case, we have from eq (5), 
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E^)=^ ik fi^iz e-^e-w-Udz, 



yielding the second term in the series for R: 

B 2 =\ik [ \ X ^) . e~ ikz dz\ik ( \ X .® e-ne-w-idt. 
With no collisions, we may write this development of R in the form 

S^Sx+^-i^Pxe- 2 ^ (12) 

If we place 

v=iE / /n 1 kE ) 

then differentiation and the use of eq (1) show immediately that 



where v=l when z^>b and 



CLV %K / o 9 9\ 

-7-= {n l —n[v l ), 

dz n x 



~7h(e' ikz +Re ikz ) 



when z<^a. This is the equation derived by Wait [1962] and which was solved by successive 
approximations. The corresponding integral equation would be 

v(fi = l+£(\n*-r&)dz, (13) 

no doubt simpler than our eq (5), but certainly less comprehensive in physical content. More- 
over, our eq (6) yields R directly as a series expansion, but Wait has shown how eq (13) may 
be solved approximately, yielding R as a series expansion only indirectly through the expansion 
of v. 

7. Examples 

As the first example, we may consider the model for which free space exists above and 
below the slab with bounding surfaces z=0 and z=h; homogeneous medium exists in the slab 
such that n 2 =l— X. It may easily be shown that the reflection coefficient is given by 

(n 2 -l)(l-e 2iknh ) 



(n+l) 2 e 2iknh -(n-l) 2 



In order to apply the second-order reflection formula (12), X is assumed to be small within 
the slab. Expanding R to the second order in X, we easily obtain 

R=iX(l-e- 2ikh )+iX 2 (l-e- 2m -2ikher 2m ). (14) 

We may also use eq (12). Since X is constant throughout the range of integration, we 
obtain 

R=^ ikX r e - 2ikz dz-~ k 2 X 2 Cdz [ e ~ 2ikz Cdt+ r e - 2ikt dt\ 

These integrals are trivial to evaluate, the result without any approximation being equal to (14). 
Secondly, we may consider a symmetrical Epstein [1930] profile, tending to free-space 
conditions above and below the layer. Let 



X- 



ae i 



0* 



(e 8z +l) 2 
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The work of Epstein (or of Budden [1961]) shows that the reflection coefficient has the value 

T(-2ik/0)T(N+2ik/0)T(l-N+2ik/l3) 

where N is given by N(N— l)+& 2 a//3 2 =0. Since a is assumed to be small, we take N=k 2 a/fi 2 
approximately, in order to calculate the first term in the development of R in terms of a. To 
this first order, we may neglect the N in the two gamma functions occurring in the numerator. 
In the denominator we write 

r (N) r (1 -N) = 7r/sin (wN)=l/N 

since N is small. Then to the first order 

R=T (-2ik/0) r (1 +2ik/l3)N 

ttN 



sin (2iirk/P) 
iriak 2 



(15) 



(3 2 sinh (2Tk/p) 

This result should also be recovered by the direct application of formula (11). This yields 

1 r m e? z e~ 2i1cz dz 

E== 2 ika J. m w+iy' 

This may be evaluated by considering the corresponding complex contour integral taken 
around the rectangle with vertices (±L, 0), (±i, 2i7r//S), the contour enclosing one double 
pole at 2=i7r/j8. The result turns out to be identical with (15). 

8. Integral Identities 

The reflection formulae (7) and (10) may be used to obtain certain interesting infinite 
integrals, when the values of 7?, X, and E are known from other considerations. 

Our first illustration involves the exact solution of a tractable model. Consider the 
exponential profile X=e aZ jk 2 with no collision frequency. The work of Budden [1961] shows 
that 

_ r(-2tfc/oQ 

' T(2ik/a) 
E(z)=— H&,a<fte""fa) 



- iv- 1 T(2ikfa) (e '*< 2 /a) - 2tk/a 
Direct substitution into eq (10) and simplification yield 

J CO 
e**e- m H® h , a (2ie aZ/2 /a)dz. 
- CO 

Finally, let 

v=2k/a, t=2e aZ/2 /a, 

where v is real and positive and arg £=0 along the positive real axis. The integral reduces to 
the result containing the single parameter v: 



-vw- l T(-iv)e* vl2 2 l - iv = C t l ~ iv H${i 



t)dt. 

Our second illustration concerns an approximate solution of eq (1). Let n 2 vary 
monotonically from 1 when z= — oo to — °o when 2= + oo ; in such a way that E tends to the 
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form e~ ikz -{-Re ikz for large negative z. Evidently n 2 possesses a zero at some real value of z; 
we shall suppose that this transition point occurs at the origin 0. Provided that the value of 
k is large enough (that is, at the higher frequencies), many authors (see for example [Heading, 
1962, page 29]) have shown that an approximate solution for E along the whole real height 
axis is given by 

E=Cn~ 1/2 f Cmlz) m Ai |7| e^' 2 k r^YH, (16) 

where Ai denotes the standard Airy integral and Cis a constant. This result is a generalization 
of the simpler W. K.B.J, solutions 



n 1/2 exp 



( ± ik\ ndz \ 



solutions that cease to be valid near the transition point 2 = 0, and that require connection 
formulae to effect a suitable connection across z=0. Equation (16), however, suffers from 
none of these disadvantages; this solution is valid even at the transition point. 
To be specific and consistent, let 

arg 7i = for 2<0, 

arg n= — }W for 2>0, 

arg z = for 2>0, 

arg z=w for z<^0. 

When £>0, we use the standard asymptotic expression 

Ai (z)-j7r- 1/2 2- 1/4 exp (-fz 3/2 ). 

Since n 2 is negative, we must write n=e~ iT/2 m J where m>0; hence 



-| Cw- ll2 e w (J k) m m- 112 exp (-k PmdA 

an exponentially evanescent solution as 2-^>oo . This solution therefore satisfies the boundary 
condition there. 

When s<0, let z=e i7C $, where f>0, yielding 

E=Cn- 1 ' 2 (e ir f m/A 1/6 Ai [7| e^ ,2 k Vne tv d^\ 

=Cn -i/2 e iH<> f r nd r y /6 Ai r /| k p wdf Y /8 i 

We now quote the standard asymptotic expression 

Ai (-x)-7r" 1/2 x- 1/4 cos (fx 3/2 -i7r) = |7r- 1/2 ^- 1/4 [exp (iix 3/2 -iiri)+i exp (-f ix^-iri)], 
yielding 

E^Cn- 1/2 e M ( r nd{\ l '"-\ t~ 1/2 (J k f wdA * 

xlexpri&j n(-fid£—-Trij+i exp(-iki n(— f)df— 4^) 

1 /% \-l/6 

^T- 1/2 Ce f ^ l^k) 7T 1/2 [exp (. . .)+% exp (. . .)]. 
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In order to introduce the incident and reflected waves, we note that if ft is large and 
positive, then 

f w(-r)rff= | r n(-r)<*r+ f V, 

Jo Jo JTj 

since n=l below z= — f x . Hence 

Jo Jo Jo 

where the upper limit is replaced by °° since the integrand must vanish at these additional low 
heights. It follows that as n-^1, 

#~i T-wCe*''* (J kV 1 [e~ ikz exp fik f " K-f)-l]^~*-7r^ 

+ie««c^(-itj"[n(-f)-l]df-iirt)]- 

In order that the incident field should be of unit amplitude, we choose 



It follows that 



B=iexp(-2ttf [n(-{-)-l]rfr)» 

, "(r* rf2 )'" Ai [(i e, " 2i /»r 



X=l-7( 2 , 



1 T -i/i e «r/a Q^'^exp^'^-f)-!]^-!^ 

Final substitution into eq (10) yields after cross multiplication 

^"^"''"(l*)' 1 " ^P (-«J"[n(-J-)-l]«if) 

=i it f " (l-7i 2 )e- a2 7i-' /2 ( f l ndzY\i |7§ e W2 /t f '»wi«Y""| da. 

This formula embraces the buildup of the reflection coefficient from every layer of a slowly 
varying medium of the type postulated. 

9. Reflection and Transmission Formulae for Anisotropic Ionospheres 

The relevant equations leading to the analysis of this section are given, for example, by 
Budden [1961, chapter III]. The differential equation for the electric field E within an aniso- 
tropic ionosphere is 

curl curl E=PE+FME, (17) 

where the 3X3 susceptibility matrix M is given by 

S-U 2 +l 2 Y 2 inYU+lmY 2 imYU+lnY 2 ^ 



Ms 



X 

V(U 2 -Y 2 ) 



-inYU+lmY 2 -U 2 +m 2 Y 2 UYU+mnY 2 

m 

{.-imYU+lnY 2 -UYU+mnY 2 -U 2 +n 2 Y 2 J 
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Here, X, Y, Z are the usual dimensionless parameters describing the ionosphere, £7=1— iZ, 
and — (I, m, n) are the direction cosines of the earth's magnetic field. The vector ME is pro- 
portional to the electric polarization vector. 

In evaluating curl E, we choose the customary coordinate system in which b/dx=—ik 
sin 6, b/dy=0, and d/ds is replaced by a prime. In component form, eq (17) may be written as 



E'J+ik sin 6E Z = - k 2 E x -k 2 (ME) x , 

E f y '-k 2 sin 2 dE y =-k 2 E y -k 2 (ME) y , 

ik sin 6E x -k 2 sin dE z =-k 2 E z -k 2 (ME) z . m 



(18) 



We now multiply these three equations by e~ ikz cos e and integrate from a to b, namely 
over the whole height range in which ionization occurs, assuming free space for z^>b and z<^a. 
Noting that 



We' 



ifcz cos QJry TPj — ikz cos i 



e dz=Ee 



b r*b 

+ik cosd Ee~ ikz cos d dz, 

a J a 



fb \b f*b 

E"e~ ikz cos d dz=(E'+ik cos 6E)e~ ikz cos e \ —k 2 cos 2 6 Ee~ ikz cos d dz, 
we obtain from eq (18) 

(E' x +ik cos 6E x +ik sin dE z )e~ ikz cos 6 \=-k 2 f [sin 2 dE x +(ME) x -sin 6 cos dE z \e~ ikz cos d dz, 



(E' y +ik cos 6E y )e- ikz cos | & =-F f (ME)^- ftl cos e dz, 

\a J a 

b nb 

1 =— k 2 \ [cos 2 dE z -s'm 6 cos 6E X + (ME) z ]e- ikzcoad 



(19) 



ik sin dE x e- ikzcoad 



Multiplying the third of these equations by tan 6 and adding it to the first, we obtain the 
simpler equation 



(E x +ik sec dE x +ik sin QE z )e~ ikz cos 



-ifc 8 f & [(ME) a; +taii 0(ME)Je- fe ' 



•<fe. (20) 



When the incident field is horizontally polarized, let 

E x =cosd ± R\ } e ikz cosd , 

E y =e- ikz cos 9 +±B±e ikt cos e . 

E z =sind ± R\\e m cosd 
for 2<a and for z>&, 

£ , a: =cos0j.r||e- < * z cos(? , 

W . T , ,,-iJcz cos 

E z =—smdj_T {] e- ikz cos9 , 
the x- and f-f actors being suppressed. When the incident field is vertically polarized, let 

E x =cos 6e~ ikz cos *+cos d\\R\]e ikz cos e , 

E y = ll R x e m coae , 

E z =—sii\ 6e~ ikz cos *+sin 0\\R\\e ikz cos B 
for 2<a and for 2>6, 
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E x =cos6\ { T\\e- m cos e , 
Ey=\\T L e' ikzC0B \ 

E 2 =—smdiiTi\e~ ikz C08d . 

Forming an arbitrary linear combination of these two fields, we write 
E X =B cos6e- ihz cos e +(A ± R\i+B\iRii) cos 6e ikz cos e , 
E y =Ae~ ikz cos e +(A ± R ± +B {] R ± )e ikz cos 9 , 



E ? =-B 



sin Be 



— ikz cos 6 



+ (A ± Rn+B il R{) sm9e ikz coa \ 



(21) 



for s<a and for 2>6, 



(22) 



#*= (4jT| + £||Tn) cos Be ~ ikz ms % 

Er=-(AiTi+B x T$ sin *«-'*»-•. J 

Equation (19) now reduces to 

-2i4 cos 6»(^l 1 i? ± +£||ff x ) = -/t 2 1 (M.E) v e- ikieos 'dz, 

J a 

while (20) takes the form 

-2ifc(il i .J2 i +jB|B B )=-i 8 J [(ME^+tan e(ME) z ]e" ikzcoa6 dz. 

The physical interpretation of these integrals at oblique incidence is demonstrated by 
considerations given by Westcott, [1962 b, c] and in further papers to be published in this series- 

In particular, if the ionized layer is weak so that the field in the medium may be replaced 
by the incident field as the first approximation, we obtain 

A ± B JL +B\\B ± =~ ik sec d C(M 2l B cos d+M 22 A-M 22 B sin d)e~ 2ikz cos 9 dz, 

1 C b 
-AJ2|| + B|i#ii=—2 ik (M n B cos d+M 12 A-M n B sin d+M 31 B sin 6 

+M Z2 A tan 6-M zz B tan 6 sin e)e- 2ikgcoa0 dz. 



Hence, if Z is constant, we obtain 

1 C b 1 TP—m 2 Y 2 C n 

L R ± =-±iksecd\ M 22 e~ 2ikz cos 6 dz=^ flfoec d ^(U 2 -™ \ Xe ~ 



m 2 Y 2 r 

''-Y 2 ) J a 



i c b 

l|i?j_=— ik sec 6 (M 2 i cos <9-M 23 sin ff)e~ 2ikz cos °dz 

w J a 



1 ., n (inYU-lmY 2 ) cos 6>+(^FD r +?7inF 2 ) sin 6 f* _ 2ik2COsBJ 
=^sec0-^ U(U 2 --Y 2 ) ~J« 2 '** C08 ^, 



L 2?H=— ^ <£& P (M 12 +M 32 tan 0)<r 2 '* 2 cos ^ 



1 . 7 (-mYF-ZmF 2 ) +tan 9(ilYU-mnY 2 ) f» 

=2 lk — ?7(C7 2 -r 2 ) — ' 



r 



l r & 

lji?ll=— ik (M„ cos 0— M 13 sin 0+M 31 sin 0-M 33 tan (9 sin 6)e- 2ikz cosd dz 



_1 ., U 2 cos 20 sec 0-(Z 2 cos 6—n 2 tan (9 sin d)Y 2 -2im sin (#77 f & ^ 
~2 ~ U(U 2 -Y 2 ) 
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J a 



Moreover, we easily deduce the formulae for the isotropic case with F=0 at oblique incidence, 

± R ± =~ ik sec 6U~ l f Xe' 2ikz ™* e dz, 



HJ2n=4 ik cos 2d sec dU ~ 1 f -X^"""* cose dz. 



The transmission coefficients may be found by multiplying eq (18) by e ikz cos e and inte- 
grating over the height range in which ionization occurs. We then multiply the third of the 
resulting equations by tan 9 and subtract it from the first, yielding 



(E r y -ik cos 6 E y )e ikzcos 
and 

{E' x -ik sec 6E x +ik sin 6E z )e ikz cos 6 



b =-k 2 r (ME) y e- ikzc08d dz, 

a J a 

--=-¥ C [(ME^+tan 6(ME) z ]e ikz cosd dz. 

J a 



Inserting the values (21) and (22) of the field above and below the range of ionization, 
we obtain 

A ± T ± +BnT±=A-~ ik sec d f (ME)^* 2C0S * efe, 

A L T\\ + B\iTi\=B-^ ik f [(ME) x +tan d(ME) 2 ]e ilcz cos e dz. 

Finally, we deduce the values of the transmission coefficients to the first order, when the 
field throughout the medium is replaced by the incident field: 

1 r b 

±T±=1— s" ^ sec e M 22 dz, 

1 C b 

I|7 7 jl= — o ^ sec e (^21 cos 6— M 23 sin 0)^2, 

x^il=— | ii£* (M 12 +tan 0M 32 )d2, 

i r & 

II jTq=1 — — iAr (Mu cos 0— M i3 sin 0+M 3 i sin 0— M 33 tan (9 sin 6)dz. 
In the isotropic case, these formulae reduce to 

± T ± =l+liksecdU~ 1 f Xdz, 



+1 ik cos 20 sec At/" 1 f Jfrfe. 

J a 



10. Comparison of the Two Methods 

We should conclude with a note on the difference between the results of Wait [1962] and 
of the method presented in this paper. Considering the isotropic model at vertical incidence, 
the first-order reflection coefficient produced here is given by (11), namely 



\ikC Xe~ 2ikz dz (23) 



R ~2 

for zero collision frequency. On the other hand, the first approximation to the solution of 
the integral eq (13) is obtained by taking v=l, yielding 
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v(fi=l+ik C (n 2 -l)dz=l+ik C Xdz. 
Hence, when f=a, 

e- ika +Re ika ^ J a 



R=±ike~ 2 



P ' Xdz (24) 

J a 



when the left hand side is expanded by the binomial since \R\ is small. 

It is obvious that this result is only valid for a thin layer, in which case the phase factor 
e -2uz m tne integrand of (23) is replaced by e~ 2i1ca . Hence, for the homogeneous slab discussed 
in section 7, eq (24) yields (when a = 0, b = h), 



lJ 2 ik Sl 



* Xdz=^ikKX. 



The exact first term (14) reduces to this only when kh is small. On the other hand, it should 
be pointed out that Wait's expansion for the homogeneous slab is of a different kind, valid 
when |M V / (1-X)|«1. 
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